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Qh! Abstract 



We present total and differential cross sections for Wbb and Zbb production at the CERN Large 



' Hadron Collider with a center-of-mass energy of ^/s = 14 TeV, including Next-to-Leading Order 

, (NLO) QCD corrections and full bottom-quark mass effects. We also provide numerical results 

I obtained with a center-of-mass energy of ^/s = 10 TeV. We study the scale uncertainty of the 

' total cross sections due to the residual renormalization- and factorization-scale dependence of the 

o . 

. truncated perturbative series. While in the case of Zbb production the scale uncertainty of the 

total cross section is reduced by NLO QCD corrections, the Wbb production process at NLO in 
^ ' QCD still suffers from large scale uncertainties, in particular in the inclusive case. We also perform 

a detailed comparison with a calculation that considers massless bottom quarks, as implemented 
in the Monte Carlo program MCFM. The effects of a non-zero bottom-quark mass (m^,) cannot be 
neglected in phase-space regions where the relevant kinematic observable, such as the transverse 
momentum of the bottom quarks or the invariant mass of the bottom-quark pair, are of the order 
of rrife. The effects on the total production cross sections are usually smaller than the residual scale 
uncertainty at NLO in QCD. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) at CERN (Geneva, Switzerland) is scheduled to start 
operation by the end of 2009. One of the most important items on its agenda is the in- 
vestigation of the mechanism of electroweak symmetry breaking (EWSB), in particular the 
discovery of one or more Higgs bosons. Once discovered, the measurement of the Higgs 
bosons' properties will be crucial to unravel the underlying EWSB mechanism. The pro- 
duction of a weak gauge boson, (from now on indicated simply as W unless differently 
specified) or Z, with a pair of bottom (6) quarks, contributing to both the W/Z + 1 6-jet and 
W/Z + 2 6-jets signatures, represents both an interesting Standard Model (SM) signal and 
one of the most important background processes to many Higgs-boson production channels. 

The cross sections for W and Z boson production with bottom quarks has been measured 
at the Tevatron pp collider at Fermilab {^/s = 1.96 TeV) by both the CDF and 
DO [3, J] collaborations. These measurements will continue with increased precision, which 
will provide a unique opportunity to test and improve the theoretical description of heavy- 
quark jets at hadron colliders by performing a thorough comparison between the Tevatron 
experimental data and existing theoretical predictions. Studying the same cross sections 
in the very different kinematic regimes available at the LHC pp collider will then be of 
great interest and will represent a crucial test of our understanding of QCD at high-energy 
colliders. 

Moreover, the production of aW and Z boson with one or two b jets represents a crucial 
irreducible background for several Higgs-boson production channels at the LHC. Wbb pro- 
duction is an irreducible background to WH associated production followed by the decay 
H bb. This is a difficult channel but theoretically very interesting, since it can play 
a very important role in measuring the 6-quark Yukawa coupling for a light Higgs boson 
at the LHC. Analogously, for a light Higgs boson, Zbb production is a background to ZH 
associated production followed by the decay H — > bb. More importantly, for heavier Higgs 
bosons, if the b quarks in Zbb production decay leptonically, the Zbb production process 
can be a background to the inclusive production of a Higgs boson followed by the decay 
H — s> ZZ with each Z boson decaying leptonically. Finally, Zbb production is a background 
to searches for Higgs bosons with enhanced 6-quark Yukawa couplings, produced in Hbb 
associated production followed by the decay H /^"""/i^ 01 H t^t~ [Sj. 
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All Higgs-boson production channels have been calculated including at least next-to- 
leading order (NLO) QCD corrections (see, e. g., Ref. for a recent review). The hadronic 
cross sections for gg ^ H and associated WH and ZH production are also known at next- 



to-next-to- 
and Refs. 



eading order (NNLO) in QCD, Refs. 



B, y, y, Q, 



11 



12 



13 



14 



16 



li 



19l . l20|, respectively. The NLO electroweak corrections to these processes 



have been calculated as well (see Ref. 



2l| and references therein {gg H) and Ref. 22 1 



{WH,ZH)). The cross section for Hbb associated production is known at NLO in QCD 



including full 6-quark mass effects 
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25 



26|. 



The production of a. W or a. Z boson with up to two jets, one of which is a 6 jet^ 



las been 



calculated including NLO QCD corrections in the variable-flavor scheme (VFS) [27l. l28l. |29|. 
while the production of a or Z boson with two b jets has been derived at NLO in QCD 
le fixed-flavor scheme (FFS), first in the massless 6-quark approximation 



using t 
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34 



35( 1 and more recently including full 6-quark mass effects 



36 



37 



38 



31 



. In the FFS only 



massless- quark densities are considered in the initial state, hence the alternative name of 4- 
fiavor- number scheme (4FNS), while in a VFS an initial-state 6-quark density is introduced, 
hence the alternative name of 5-fiavor- number scheme (5FNS). The two schemes amount to 
a different ordering of the perturbative series for the production cross section: in the 4FNS 
the perturbative series is ordered strictly by powers of the strong coupling as, whereas in 
the 5FNS the introduction of a 6-quark parton distribution function (PDF) allows to resum 
terms of the form a" ln(m^/M^)'^ at all orders in Og (for fixed order of logarithms m), where 
M represents the upper integration limit of the 6-quark transverse momentum and can be 
thought to be of the order of Mw or Mz- While the two approaches can give very different 
results at the leading or lowest order (LO) in QCD, starting at NLO in QCD they have been 
shown to be consistent within their respective theoretical uncertainties for both H + 1 b-jet 



production 



39 



40 



4l| (for a brief review see also Ref. 26|) and single-top production 42 1. 



Recently, ly + 1 6-jet production has been calculated by consistently combining both NLO 
4FNS and 5FNS calculations 43|]. Since the relevance of the logarithms resummed in the 
VFS approach varies with the kinematic regime considered, combining the two calculations 
improves the accuracy of the theoretical prediction for + 1 6-jet production. A similar 
study is currently in progress for Z + 1 6-jet production 4^. Improving the predictions for 
Z + 1 6-jet production will be particularly relevant at the LHC, where this process allows 
for a direct determination of the 6-quark PDF, to be used in the prediction oi H + 1 6-jet 
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production, a discovery channel for beyond-the-SM Higgs bosons with enhanced 6-quark 
Yukawa couphngs. 

In this paper we provide results for Wbb and Zbb production at the LHC, keeping the 
W and Z boson on shell and with both b jets tagged in the final state, i. e. we focus on 
the W/Z + 2 6-jet case. We include NLO QCD corrections and full 6-quark mass effects. 
The corresponding results for Wbb and Zbb production at the Tevatron have been presented 



in Refs. 



36 



37, 



The details of the NLO QCD ca 



36 



culations used to compute the NLO 
. In this paper we focus on 



37 



38 



cross sections in this paper can be found in Refs. 
presenting results for the total production cross sections and kinematic distributions that are 
of interest to LHC physics. Using the MCFM package 45|], we compare our results to those 
obtained by neglecting the 6-quark mass at NLO in QCD. A non-zero 6-quark mass mainly 
affects phase space regions where the relevant kinematic observable, such as the transverse 
momentum of the bottom quarks or the invariant mass of the 6-quark pair, are of the order 
of rrih. Indeed, this is the reason why 6-quark mass effects cannot be neglected in studies 
of W/Z + ra-jet production (n = 1,2) with at least one b jet, as has been discussed for the 



n 



2 case in Refs. 



28 



291] and for the n = 1 case in 431 ]. 



The paper is organized as follows: in Section [TTl we briefly describe our choice of input 
parameters, cuts, jet-identification algorithm, and observables, while we present numerical 
results and a discussion of the 6-quark mass effects for Wbb and Zbb production in Section [XTTl 
and Section IIVI respectively. Section |V] contains our conclusions. 



II. GENERAL SETUP 

The results for both Wbb and Zbb production presented in this paper have been obtained 
for the LHC pp collider running at a center-of-mass energy of either ^/s = 10 TeV or ^/s = 
14 TeV. While we only provide kinematic distributions obtained with a center-of-mass energy 
of ^/s = 14 TeV, we compare the total cross sections obtained for center-of-mass energies of 
^/s = 10 TeV and ^/s = 14 TeV in Tables HI and HTl The mass of the bottom quark is taken to 
be rrify = 4.62 GeV. Results in the massless 6-quark approximation have been obtained using 



the MCFM code (version 5.4) 45|] . The top-quark mass, entering in the virtual corrections, is 
set to nit = 172.6 GeV. In the case of Wbb production we use = 80.44 GeV, while for Zbb 
production we use Mz = 91.1876 GeV and derive Mw from the relation Mw = Mz cos 9w 
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In both cases we assume sin^ 9w = 0.223. We work in the electroweak input scheme and 
replace the fine structure constant a(0) = e^/(47r) by a{G^) = fj_M^ sin^ 9^ with the 
Fermi constant = 1.16639-10-5 GeV-^ The W-h oson couphng to quarks is proportional 
to the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. We use non-zero CKM matrix 
elements for the first two quark generations, Vud = Vcs = 0.974 and Vus = Vcd = 0.227, while 
we neglect the contribution of the third generation, since it is suppressed either by the 
initial-state quark PDFs or by the corresponding CKM matrix elements. 

The LO results are based on the one-loop evolution of and the CTEQ6L1 set of 
PDFs 46(1, with a^^{Mz) = 0.130, while the NLO results use the two-loop evolution of 
and the CTEQ6M set of PDFs, with a^^^{Mz) = 0.118. In the calculation of the parton 
luminosity we assume five light fiavors in the initial state, but we have verified that including 



the fe-quark PDF has a neg. 



igible effect (< 0.1%) on the W/Zbb cross section. We implement 



the kx jet algorithm [47|, |48|, |49|, l50| with a pseudo-cone size of i? = 0.7 and we recombine the 



parton momenta within a jet using the so called covariant E'-scheme |48[. We checked that 
our implementation of the jet algorithm coincides with the one in MCFM. We require all 
events to have a 6-jet pair in the final state, with a transverse momentum larger than either 
15 GeV or 25 GeV, in order to study the dependence on the 6-jet transverse-momentum 
cut. We also require that the pseudorapidities of both b jets satisfy \ri^'^\ < 2.5. We impose 
the same pt and |?7| cuts also on the extra jet that may arise due to hard non-coUinear 
real emission of a parton, i. e. in the processes W/Zbb + g oi W/Zbb + q{q). This hard 
non-collinear extra parton is treated either inclusively or exclusively. In the inclusive case 
we include both two- and three-jet events, while in the exclusive case we require exactly two 
jets in the event. Two-jet events consist of a 6-jet pair that may also include a final-state 
light parton (gluon or quark) due to the applied recombination procedure. On the other 
hand, three-jet events consist of events containing a 6-jet pair plus an extra light jet. We 
notice that, at NLO in QCD, all jets in three-jet events consist of a single parton. 

For both Wbb and Zbb production we provide results for the total cross section (cr) 
and the following kinematic distributions: da/dp^\ da/dpjf'', da/dp'!^^^ , da/dri^'\ da / dt]''''^'' , 
da/drf^/^ , da/dm^^l, and da/dR^'^, where p^', p^^'', and p^^^ are the transverse momenta 
of the leading b jet (i. e. leading in pt), the subleading b jet and of the W or Z boson, 
^b,i^ ^b,si^ ^w/z ^Yi^Q corresponding pseudorapidities, m^g is the invariant mass of the 
bb pair, and Rjji is their relative separation in the pseudorapidity-azimuthal angle plane. 



5 



70 

60 

50 

^40 

D 30 

20 

10 


50 

40 



X) 30 



20 



10 



0.5 



p^>15GeV |Ti|<2.5 





H h 



5 — 



p^>25GeV |Ti|<2.5 



W bb 



0.5 



W bb 



1 ^ ^ n 

X — »< NLOInc 
V — V NLOExc 
o — c LO 




Total 




FIG. 1: Dependence of the LO (black, solid), NLO exclusive (blue, solid), and NLO inclusive (red, 
solid) total cross sections for W~^hh and W~bb production on the renormalization/factorization 
scales, including full 6-quark mass effects, when fi = fir = f^f is varied between fio/2 and 2fj,Q (with 
^0 = M\Y + 2mh). We also show the individual channels, qq' (dash-dotted) and qg + qg (dotted), 
for the inclusive (red) and exclusive (blue) cases. 



Rf^l = ^(ryb,' - r/M)2 + (0fe,« _ 0fe,.i)2, 

III. Wbb PRODUCTION 

At tree level, the production oi aW boson with a pair of bottom quarks consists of just 
one process, qq' Wbb. In order to compute this process at NLO in QCD one needs to 
include one-loop virtual corrections to qq' — > Wbb as well as all real radiation corrections 
with up to one extra parton in the final state, i. e. qq' — » Wbb + g and qg{qg) Wbb + q'{q'). 



Details of the calculation have been given in Refs. 36|, |38|| and will not be repeated here. 



We note that, contrary to a pp collider like the Tevatron, at a pp collider like the LHC 
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TABLE I: LO, NLO inclusive, and NLO exclusive cross sections for W^bb and W~bb production 
at both y/s = 10 TeV and -y/s = 14 TeV, for two different values of the b-jet transverse-momentum 
selection cut, and for a zero (obtained with MCFM) and non-zero 6-quark mass. The central 
values correspond to Hr = /.tj = /^o = Mw + 2mfe, while the upper and lower bounds represent the 
maximal upper and lower variation obtained when varying fir = l^f between /io/2 and 2fXQ. 
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the cross sections for W^hb and W~hh are different, because the two processes depend on 
different quark/antiquark PDFs that are not symmetrically distributed between the two 
incoming nucleons. In the following we will provide results separately for both production 
processes. 

In Fig. [T] we illustrate the renormalization- and factorization-scale dependence of the 
LO and NLO total cross sections obtained for a massive h quark, when /i = /i^ = /// is 
varied between /io/2 and 2/io, with /iq = M\y + 2mh. We immediately notice that the 
impact of NLO QCD corrections is very large, in particular in the inclusive case, where they 
increase the LO cross section by a factor between two and three depending on the scale. 
We also notice that the scale dependence of the NLO cross section is worse than {inclusive 



case) or comparable to [exclusive case) the scale dependence o: 



is different from what has been observed for the Tevatron 
out in a calculation with massless bottom quarks 



36 



the LO cross section. This 
381 ] ■ and was first pointed 



35| . It is just a reminder of the fact 
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FIG. 2: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
transverse momentum distributions for the 6 jet with the leading (left hand side) and subleading 
(right hand side) transverse momentum in W'^bb (upper plots) and W~bb (lower plots) production. 
The lower window shows a bin-by-bin K factor, for the inclusive (red, solid) and exclusive (blue, 
dot-dashed) cases. 



that, at a given perturbative order, the uncertainty due to the residual renormalization- and 
factorization-scale dependence may underestimate the theoretical uncertainty due to missing 
higher-order corrections. A realistic determination of this uncertainty is usually much more 
complex and requires a thorough understanding of the perturbative structure of the cross 
section, in particular at the lowest orders of the perturbative expansion. In Wbb production 
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FIG. 3: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
transverse momentum distributions for the (left hand side) and W~ (right hand side) bosons 
in W'^bb and W~bb production respectively. The lower window shows a bin-by-bin K factor, for 
the inclusive (red, solid) and exclusive (blue, dot-dashed) cases. 

the NLO QCD corrections introduce a new, numerically important production channel not 
present at LO, as will be discussed in more detail below. Therefore, only at NLO the scale 
dependence of the cross sections starts to be a meaningful measure of the behavior of the 
perturbative expansion. 

In order to understand better the behavior of the W~^bb and W~bb cross section, we 
also show in Fig. [1] the scale dependence of the individual parton-level channels. We notice 
that although the qg{qg) W^bb + q'{q') channel appears for the first time at NLO in the 
perturbative expansion of the Wbb cross section, it is actually a tree level contribution and, 
as such, introduces a large scale dependence in the calculation that will be moderated only 
by adding (still unknown) NNLO corrections. The reason why this becomes so evident at 
the LHC, while it is not at the Tevatron, is because the qg{qg) W^bb + q'{q') channel is 
enhanced by the correspondingly large initial-state gluon PDF. The NLO Wbb total cross 
section is particularly affected by this process because there is no gluon-initiated process at 
LO. Finally, the impact of the qg{qg) W^bb+q'{q') channel on the scale dependence of the 
total cross section is larger in the inclusive than in the exclusive case because the exclusive 
cross section by definition discriminates against processes with more than two jets in the 
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FIG. 4: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
pseudorapidity distributions for the b jet with the leading (left hand side) and subleading (right 
hand side) transverse momentum in W'^bb (upper plots) and W~bb (lower plots) production. The 
lower window shows a bin-by-bin K factor, for the inclusive (red, solid) and exclusive (blue, dot- 
dashed) cases. 



final state. Figure [U also shows the effect of lowering the cut on the transverse momentum 
of the b jets. Lowering the cut from > 25 GeV to > 15 GeV almost doubles the cross 
section and can therefore be a crucial factor in deciding which selection cuts to use for b 
jets. 

In Table [T] we present both LO and NLO total cross sections separately for W~^bb 
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FIG. 5: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
pseudorapidity distributions for the (left hand side) and W~ (right hand side) bosons in 

W~^bb and W~bb production respectively. The lower window shows a bin-by-bin K factor, for the 
inclusive (red, solid) and exclusive (blue, dot-dashed) cases. 



and W~bb production, including our estimate of the scale uncertainty due to the resid- 
ual renormalization- and factorization-scale dependence. We provide results for both 
> 15 GeV and > 25 GeV, and for both center-of-mass energies, ^/s = 14 TeV 
and y/s = 10 TeV. We also include the corresponding set of results obtained with MCFM 
for mi, = 0. Comparing the results from the massless approximation and our results with full 
6-quark mass dependence, one observes that for p'^ > 15 GeV the massless approximation 
overestimates the total cross section by about 10%, while for p^ > 25 GeV the difference is 
a milder 3% (both for ^/s = 10 TeV and ^/s = 14 TeV). As expected, the more inclusive 
the treatment of the b jets the more important the 6-quark mass effects become. This, for 
example, explains why for a complete NLO treatment of W + 1 6-jet production, the con- 
tributions from the Wbb production process must be calculated using the full 6-quark mass 
dependence, as discussed in [43|. 

In Figs. (I2j)-(l5]) we show the transverse momentum (pt) and pseudorapidity (77) distri- 
butions for the b jet with the leading and subleading transverse momentum and for the W 
boson for both W^bb and W~bb production. The results for these distributions, as well as 
for all other distributions presented in this paper, have been obtained with y/s = 14 TeV 
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m,, [GeV] 



FIG. 6: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) bb- 
pair invariant mass distributions in W~^bb (left hand side) and W~bb (right hand side) production. 
The lower window shows a bin-by-bin K factor, for the inclusive (red, solid) and exclusive (blue, 
dot-dashed) cases. 



and by assuming > 15 GeV and \ri''\ < 2.5. The upper panels of each figure show the 
LO, NLO inclusive and NLO exclusive distributions, while the lower panels show the ratios 
da^^lo/do'LO and da'^lQ/dcLo, thereby providing a bin-by-bin K factor. In each figure we 
show the results obtained with the central scale choice fir = fJ^f = fio = + 2mb. The 
distributions for the invariant mass of the two h jets (m^s) and for their relative separation in 
the pseudorapidity-azimuthal angle plane (-Rbg) are shown in Fig. ([6]) and ([7]), respectively. 

Clearly, the NLO QCD corrections largely affect the kinematic distributions, resulting in 
considerable changes in their shapes, both in the inclusive and exclusive case. Figure ([2]) 
shows that at NLO in QCD the production of h jets at large p\ is consistently suppressed 
in the exclusive case, while in the inclusive case the cross section for the production of a 
leading h jet is enhanced, yielding K factors of about 2 in the low region and of about 3.5 
at large p^. Similar features can be observed in the distributions of Fig. ([3]). We note 
that this may have an impact especially on assessing the effects of experimental triggers on 
the lepton coming from the H^-boson decay. The pseudorapidity distributions of the leading 
and sub-leading h jet of Figs, (jlj) and ^ in the exclusive case are enhanced by the NLO 
QCD corrections but their shape is barely affected, while in the inclusive case the increase 
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'^bb ^b 

FIG. 7: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) Rj^i 
distributions in W'^bb (left hand side) and W~bb (right hand side) production. The lower window 
shows a bin-by-bin K factor, for the inclusive (red, solid) and exclusive (blue, dot-dashed) cases. 

of the rj^ distribution of the leading b jet is more pronounced in the central region. The 
increase of the pseudorapidity distributions of the W bosons at NLO QCD, on the other 
hand, is more pronounced in the forward regions in both the inclusive and exclusive cases. 
Finally, we point out that the large positive corrections to the Ryi, distributions of Fig. ([7]) 
in regions with low and large values of Rf,i are especially pronounced in the inclusive case. 
Although these changes in shape are commonly seen in NLO QCD computations involving 
two or more jets, the large effects observed here in the inclusive specific feature of 

Wbb production, originating from the tree-like g(7-initiated contribution to the NLO QCD 
corrections, as one can easily deduce from comparing with the exclusive K factors. 

Finally, we compared the distributions presented here to the ones produced by MCFM 
in the massless 6-quark approximation. We agree very well, with only small (of the order 
of the change in the total cross sections shown in Table [I]) but noticeable deviations in 
regions where relevant kinematic observables become small, i. e. comparable to nih. As an 
example, we show in Fig. [8] the comparison of the LO and NLO m^^g and Ri,i distributions for 
W~bb production obtained from the massive and massless 6-quark calculations. Most of the 
difference between the massless and massive 6-quark cross sections comes from the region of 
low invariant mass m^g, and are more pronounced for small values of Rf,i, both at LO and 
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FIG. 8: LO (black), NLO inclusive (red) and NLO exclusive (blue) m^j (left hand side) and Rf^i 
(right hand side) distributions for W^bb production derived from our calculation with rrih ^ 
(LO: dashed, NLO inclusive: solid, NLO exclusive: dash-dotted) and from MCFM with nib = 
(LO: double-dashed/dotted, NLO inclusive: dashed/double-dotted, NLO exclusive: dotted). The 
lower window shows the ratio of the distributions for massive and massless b quarks, da{'mij ^ 
0)/da{mi, = 0) (LO: dashed, NLO inclusive: solid, NLO exclusive: dash-dotted). 



at NLO, where the cross sections for rrib ^ are consistently smaller than the ones with 
= 0. This may indicate that a resummation of large logarithmic corrections may be in 
order when the two b jets become collinear. As can be seen by comparing the ratios of the 
LO and NLO cross sections in Fig. [HI the impact of a non-zero 6-quark mass is almost not 
affected by including NLO QCD corrections and can be taken into account by rescaling the 
NLO result for massless bottom quarks with the ratio of the LO cross sections as discussed 
in [36 1 . 



IV. Zbb PRODUCTION 

At tree level, the production of a Z boson with a pair of bottom quarks consists of two 
channels, namely qq Zbb and gg Zbb. At NLO in QCD one needs to include the one- 
loop virtvial corrections to both tree-level processes as well as the real radiation corrections 



14 



Pj>15GeV lTil<2.5 



p^>25GeV lTil<2.5 



200 



150 



X> 100 



50 



1 ' 1 ' 


1 




Lob 




^-"^ 










1 , r 


.■ 1 


0.5 1 


2 







80 



60 



40 



20 — 



NLO Inc 
NLO Exc 
LO 




Total 



qq 

..... qg 



0.5 



FIG. 9: Dependence of the LO (black, solid), NLO exclusive (blue, solid), and NLO inclusive (red, 
solid) total cross sections for Zhh production on the renormalization/factorization scales, including 
full 6-quark mass effects, when /i = /i^ = ^/ is varied between /io/2 and 2^q (with = Mz + 2mi,). 
We also show the individual channels, qq' (dash-dotted), qg + qg (dotted) and gg (dashed), for the 
inclusive (red) and exclusive (blue) cases. 



with up to one extra parton in the final state, i.e. qq Zhh + g, gg ^ Zhh + g, and 



Qgi^g) ~^ Zhh + q{q). Details of the calculation have been given in Refs. 37|, |38| and will 
not be repeated here. 

As done in the case of W^hb and W~hh production in Section IIIIl we start our discus- 
sion of Zhh cross sections at NLO QCD at the LHC by studying the renormalization- and 
factorization-scale dependence of the LO and NLO total cross sections. In Fig. [9] we show 
the scale dependence of the individual parton-level channels, as well as their sum, at LO and 
NLO in QCD for both the inclusive and exclusive case. As for Whb production, we notice 
the pronounced scale dependence of the NLO total cross section to qg{qg) — * Zbb + q{q). 
However, since in Zhh production both the LO and NLO cross sections also consist of a 
g g-'miiiated subprocess, the qg{qg) —>■ Zhh + q{q) channel is not dominant at NLO and its 
effect is therefore less pronounced. Indeed, the scale dependence of the exclusive cross sec- 
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TABLE II: LO, NLO inclusive, and NLO exclusive cross sections for Zbb at both ^/s = 10 TeV and 
y/s = 14 TeV, for two different values of the 6-jet transverse-momentum selection cut, and for zero 
and non-zero 6-quark mass. The central values correspond to fir = f^f = f^o = Mz + 2mh, while the 
upper and lower bounds represent the maximal upper and lower variations obtained when varying 
/ir = M/ between fio/2 and 2fiQ. 
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tion actually greatly improves at NLO in QCD, while the scale dependence of the inclusive 
one is only mildly better than at LO, but not worse as it is the case in Wbb production. 

In Table Hn we provide both LO and NLO total cross sections for Zbb production, including 
our estimate of the residual uncertainty due to only the renormalization- and factorization- 
scale dependence. We give results for two choices of cuts, > 15 GeV and p^ > 25 GeV, 
and for both center-of-mass energies, ^/s = 14 TeV and ^/s = 10 TeV. We compare the case 
of a massless and massive bottom quark and we observe that the effects of a non-zero b- 
quark mass are mild for Zbb production. Indeed, for p^ > 15 GeV the massless approximation 
mildly overestimates the total cross section by about 2 — 3%, which in the inclusive case is 
considerably smaller than the scale uncertainty. For p'l > 25 GeV the difference is basically 
gone. We remind, however, that this is not the case in more inclusive studies, like Z + 1 6-jet 
production 44 1, where, as seen in ly + 1 6-jet production {43!, one needs to consistently 
include full 6-quark mass effects. 

In Figs. (fT0|) - (fT3|) we show the transverse momentum (pt) and pseudorapidity (77) distri- 
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FIG. 10: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
transverse momentum distributions for the 6 jet with the leading (left hand side) and subleading 
(right hand side) transverse momentum in Zbb production. The lower window shows the bin-by-bin 
K factor, for the inclusive (red, solid) and exclusive (blue, dot-dashed) cases. 



butions for the b jet with the leading and subleading transverse momentum, and for the Z 
boson. As in the Wbb case, the upper parts of each figure show the LO, NLO inclusive and 
NLO exclusive distributions, while in the lower parts we provide the bin-by-bin K factors. 
We show results obtained by using the central scale choice /i^ = /x/ = /io = Mz + 2mb. As 
expected from the study of the scale dependence shown in Fig. ([9]), in the exclusive case 
the NLO QCD corrections reduce the cross sections as compared to the LO result, while 
in the inclusive case the NLO QCD corrections increase them. Figures fllUj) and f ill I) also 
show, in particular in the bin-by-bin K factors, that the NLO QCD corrections affect the 
shape of the exclusive transverse momentum distributions, with large negative corrections 
at large px where the bin- by-bin K factors become smaller than one, while the changes in 
the shape are much less pronounced in the inclusive case. Figures f|T2|) and f|T3|) show that 
NLO QCD corrections have hardly any effects on the shape of the pseudorapidity distribu- 
tions. The distributions for the invariant mass of the two b jets (m^j) and for their relative 
distance {Ri,i) are shown in Figs. (HM and (fT^ respectively. We notice that, although in a 
much milder fashion than for Wbb production, in the inclusive case NLO QCD corrections 
enhance the regions at small and large values of Rf,i. 
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FIG. 11: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
transverse momentum distributions for the Z boson in Zbb. The lower window shows the bin-by-bin 
K factor, for the inclusive (red, solid) and exclusive (blue, dot-dashed) cases. 



Finally, in Fig. [16] we compare the LO and NLO rrif^i and Rf,i distributions for Zbb 
production obtained from the massive and massless 6-quark calculations. The results with 

= have been obtained using MCFM. Most of the difference is seen in the region of low 
invariant mass mf,i, and are more noticeable for small values of Ri,i. In all regions though the 
impact of a non-zero 6-quark mass is almost not affected by including NLO QCD corrections, 
and can be taken into account by rescaling the NLO result for massless bottom quarks with 
the ratio of the LO cross sections as discussed in 371]. 



V. CONCLUSIONS 



A reliable theoretical prediction for W and Z production with b jets will be crucial for 
many Higgs-boson searches and studies at the LHC. In this paper we have presented a 
study of W and Z boson production with two b jets at the LHC including full 6-quark 
mass effects, based on the NLO QCD calculation of Wbb and Zbb production presented in 



Refs. 



36 



37 



We have presented numerical results for the total Wbb and Zbb production 
cross sections, as well as for a number of kinematic distributions interesting to LHC physics, 
for both massive (our calculation) and massless bottom quarks (as implemented in MCFM). 



18 




FIG. 12: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
pseudorapidity distributions for the b jet with the leading (left hand side) and subleading (right 
hand side) transverse momentum in Zbb production. The lower window shows the bin-by-bin K 
factor, for the inclusive (red, solid) and exclusive (blue, dot-dashed) cases. 



The Wbb cross sections at NLO QCD still suffer from a large theoretical uncertainty due 
to the unphysical renormalization- and factorization-scale dependence, which is particularly 
pronounced in the inclusive case, and further theoretical improvements are warranted. In 
the case of Zbb production, the NLO QCD cross sections are well behaved, i. e. exhibit 
only a mild residual scale dependence, in particular in the exclusive case. The shape of 
distributions is changed significantly by NLO QCD corrections, such that they can not be 
correctly described by global K-factor rescalings. This is more dramatic in the case for 
Wbb production, in view of which including resummation effects, as well as consistent NLO 
showering, might be of considerable importance. 

The fe-quark mass effects can impact the shape of the kinematic distributions, as shown 
on the example of the and -Rfeg distributions, in particular in phase space regions where 
the relevant kinematic observable is of the order of m^. Apart from these regions, however, 
these effects can be approximated by rescaling the NLO cross section for rrih = with the 
ratio o : 
in [3f 

the effect is more pronounced the smaller the applied cut. However, these effects are in 



LO cross sections for massive and massless bottom quarks as discussed in detail 
□ 

37| . The total production cross sections are reduced by 6-quark mass effects, and 
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FIG. 13: LO (black, dashed), NLO inclusive (red, solid) and NLO exclusive (blue, dot-dashed) 
pseudorapidity distributions for the Z boson in Zbb production. The lower window shows the 
bin-by-bin K factor, for the inclusive (red, solid) and exclusive (blue, dot-dashed) cases. 

most cases smaller than the residual scale dependence at NLO in QCD, especially in Wbb 
production for the inclusive case. 
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